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Using the nuclear symmetry energy that has been recently constrained by the isospin 
diffusion data in intermediate-energy heavy ion collisions, we have studied the transition 
density and pressure at the inner edge of neutron star crusts, and they are found to be 
0.040 fm~ 3 < p t < 0.065 fm~ 3 and 0.01 McV/fm 3 < Pt < 0.26 MeV/fm 3 , respectively, 
in both the dynamical and thermodynamical approaches. We have also found that the 
widely used parabolic approximation to the equation of state of asymmetric nuclear 
matter gives significantly higher values of core-crust transition density and pressure, 
especially for stiff symmetry energies. With these newly determined transition density 
and pressure, we have obtained an improved relation between the mass and radius of 
neutron stars. 



1. Introduction 

Studying the properties of neutron stars, which are among the most mysterious 
objects in the universe, allows us to test our knowledge of matter under extreme 
conditions. Theoretical studies have shown that neutron stars are expected to have 
a liquid core surrounded by a solid inner crust which extends outward to the 
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neutron drip-out region. While the neutron drip-out density p out is relatively well 
determined to be about 4 x 10 11 g/cm 3 or 0.00024 fm -3 1^1, the transition density 
p t at the inner edge is still largely uncertain mainly because of our very limited 
knowledge on the nuclear equation of state (EOS), especially the density depen- 
dence of the symmetry energy (E sym (p)) of neutron- rich nuclear matter 12141. These 
uncertainties have hampered our understanding of man y important properties of 
neutron stars I2|4|S] anc j re j a ted astrophysical observations | 3 | 4 | 6 | 7 | 8 | 9 | 10 | ll | 12 | 13 jl4| 
Recently, significant progress has been made in constraining the EOS of neutron-rich 
nuclear matter using terrestrial laboratory experiments (See Ref. f or the mo st 
recent review). In particular, the analysis of the isospin-diffusion data^ ^^^ in 
heavy- ion collisions has constrained tightly the E sym (p) in exactly the same sub- 
saturation density region around the expected inner edge of neutron star crust. The 
extracted slope parameter L — 3po dE ^^ \ p=Po in the density dependence of the 
nuclear symmetry energy was found to be 86 ± 25 MeV, which has further been 
confirmed by a more recent analysis With this constrained nuclear symmetry 
energy, we have obtained an improved determination of the values for the transi- 
tion density and pressure at the inner edge of neutron star crusts. This has led us 
to obtain significantly different values for the radius of the Vela pulsar from those 
estimated previously. Also, we have found that the widely used parabolic approxi- 
mation (PA) to the EOS of asymmetric nuclear matter gives much larger values for 
the transition density and pressure, especially for stiff symmetry energies. 

2. Dynamical and Thermodynamical approaches to the stability of 
npe matter 

The transition densi ty at the inner e dge of neutron star crusts can be det ermined 
using the dynamical | 6 | 7 | 8 | 9 | 20 | 21 | 22 | anc j thermodynamical approaches | 4 | 23 | 24 | 
as well as the approach based on the random phase approximation (RPA) | 25 | 26 | 
Our study was based on the dynamical and thermodynamical approaches. 

In the dynamical approach, the stability condition for a homogeneous n pe matter 
against small periodic density perturbations can be well approximated by| 6 | 7 | 8 | 9 121| 

Aire 2 

V dyn (k) = V + (5k 2 + > 0, (1) 

where k is the wavevector of the spatially periodic density perturbations and 
_ dfi p (d^ n /dp p ) 2 2 47re 2 

^0 — — ^ — ) Krpp 



dp p dpn/Opn ' dpe/dpe' 

l3 = D pp + 2D np C + D nn C, 2 , C 



2 t _ dp, n /dp p 
dp, n /dp n 



with fii being the chemical potential of particle type i. The three terms in Eq. {1} 
represent, respectively, the contributions from the bulk nuclear matter, the density- 
gradient (surface) terms, and the Coulomb interaction. For the coefficients of 
density-gradient terms, we use the empirical values of D pp = D nn = D np = 132 
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McV-f m 5 , wh ich are consistent with those from the Skyrme-Hartree-Fock calcu- 
lations ISlIlS. At fc mi „ = [(^) 1/2 - ^ F ] 1/2 , V dyn (k) has the minimal value of 

VdynOmm) = V + 2(4 7 re 2 /3) 1 / 2 - (3kl F I 6 ' 7 ' 8 ' 9 ' 21 ' and p t is then determined from 

Vdyn(fcmin) = 0. 

In the t hermodynamical approach, the stability conditions for the npe matter 
4 23H] 



are 

In the above, P = Pb + P e is the total pressure of the npe matter with Pb and P e 
being the contributions from baryons and electrons, respectively; v and q c are the 
volume and charge per baryon number; and p = p n — p, p is the chemical potential of 
the npe matter. It can be shown that the first inequality in Eq. ((2|) is equivalent to 
requiring a positive bulk term Vq in Eq. while the second inequality in Eq. 

is almost always satisfied. Therefore, the thermodynamical stability conditions are 
simply the limit of the dynamical one for k — > when the surface terms and the 
Coulomb interaction are neglected. 



3. Results for the transition density and pressure 

We have used in our study a momentum-dependent MDI interaction that is based 
on the modified finite-range Gogny effective interactional. This interaction, which 
has been extensively studied in our previous work E3, is ex actly t he one used in 
analyzing the isospin diffusion data from heavy- ion reactions l l^ l l^ l. 

Shown in Fig. Q] is the transition density p t as a function of the slope param- 
eter L of the symmetry energy from the MDI interaction. For comparisons, we 
have included results from both the dynamical and thermodynamical approaches 
with the full EOS and its parabolic approximation, i.e., E(p,S) = E(p,5 = 
0) + E sym (p)5 2 + 0(<5 4 ), where S = (p n — p P )/p is the isospin asymmetry, from 
the same MDI interaction. With the full MDI EOS, pt is seen to decrease almost 
linearly with increasing L in both approaches consistent with the results obtained 
in RPA 1251. Both dynamical and thermodynamical approaches give very similar 
results with the former having slightly smaller p t than the later (the difference is 
actually less than 0.01 fm -3 ), and this is due to the fact that the density-gradient 
and Coulomb terms in the dynamical approach make the system more stable and 
thus lower the transition density. The small difference between the two approaches 
implies that the effects of density-gradient terms and Coulomb term are, however, 
unimportant in determining p t . On the other hand, significantly larger transition 
densities are obtained in the parabolic approximation, including the predictions by 
Kubis using the MDI EOS in the thermodynamical approaches], especially for stiffer 
symmetry energies (larger L values). The large error introduced by the parabolic 
approach is understandable since the /3-stable npe matter is usually highly neutron- 
rich and the contribution from the higher-order terms in 5 is thus appreciable. This 



November 5, 2009 18:37 WSPC/INSTRUCTION FILE NScrust 



4 Jun Xu et al. 




L (MeV) 

Fig. 1. (Color online) The transition density pt as a function of the slope parameter L of 
the nuclear symmetry energy from the dynamical and thermodynamical approache s w ith 
and without the parabolic approximation in the MDI interaction. Taken from Ref. ^3. 



is especially the case for the stiffer symmetry energy which generally leads to a more 
neutron-rich npe matter at sub-saturation densities. Furthermore, because of the en- 
ergy curvatures involved in the stability conditions, larger factors are multiplied to 
the contributions from higher-order terms in the EOS than tha t multiplied to the 
quadratic term. These features agree with the early finding ^Hl that the transition 
density p t is very sensitive to the fine details of the nuclear EOS. According to 
results from the more complete and realistic dynamical approach, the constrained 
L limits the transition density to 0.040 fm~ 3 < p t < 0.065 fm -3 as shown in Fig. [TJ 
The transition pressure Pt at the inner edge of the neutron star crust is also 
an important quant ity that might be measurable indirectly from observations of 
pulsar glitches EE] Shown in Fig. [2] is the Pt as a function of L and p t by using 
the dynamical approach with both the full MDI EOS and its PA. Again, it is seen 
that the PA leads to huge errors for large (small) L (p t ) values. For the full MDI 
EOS, the P t decreases (increases) with increasing L (p t ) while it displays a complex 
relation with L or p t in PA. From the constrained L values, the value of Pt is limited 
between 0.01 MeV/fm 3 and 0.26 MeV/fm 3 . 



4. Improved constraint on the radius-mass relation of neutron 
stars 

The constrained values of p t and Pt have important implications in many properties 
of neutron stars I 4 | 8 | 12 | 2l ] ^ g an exam pi ej we have examined their effect on con- 
straining the mass-radius (M-R) correlation of neutron stars. The crustal fraction 
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Fig. 2. (Color online) The transition pressure Pt as a functionof L and pt by using the 
dynamical approach with and without the parabolic approximation in the MDI interaction. 
Taken from Ref . ^3. 



of the total moment of inertia AI 1 1 of a neutron star can be well approximated 

b .. i5Hirn 



AI 

T 



28irP t R 3 (1 - 1.67£ - 0.6£ 2 ) 



3Mc 2 



i 



2P t (l + 5£-14£ 2 ) 



PtrribC 2 ^ 



(3) 



where mi, is the baryon mass and £ = GM / Rc 2 with G being the gravitational 
constant. As stressed in Ref. AI / 1 depends sensitively on the symmetry energy 
at subsaturation densities through P t and p t , but there is no explicit dependence on 
the higher-density EOS. So far, the only known limit of AI /I > 0.014 was extracted 
from studying the glitches of Vela pulsar HH. This together with the upper bounds 
on P t and pt (pt — 0.065 fm -3 and P t = 0.26 MeV/fm 3 ) sets approximately a 
minimum radius of R > 4.7 + 4.0M/M Q km for the Vela pulsar. The radius of Vela 
pulsar is predicted to exceed 10.5 km should it have a mass of 1.4M Q . We notice 
that a constraint of R > 3.6 + 3.9M/M km for this pulsar has previously been 
derived in Ref. by using p t = 0.075 fm~ 3 and P t — 0.65 MeV/fm 3 . However, 
the constraint obtained in our study using for the first time data from both the 
terrestrial laboratory experiments and astrophysical observations is more stringent. 

The above constraints are shown in Fig. [3] together with the M— R relation ob- 
tained by solving the Tolman-Oppenheimer-Volkoff (TOV) equation. In the latter, 
we have used the well-known BPS EOS^for the outer crust. In the inner crust with 
Pout < P < Pt, the EOS is largely uncertain and following Ref.^Sl, we use an EOS of 
the form P = a + be 4 / 3 with the constants a and b determined by the total pressure 
P and total energy density e at p out and pt- The full MDI EOS and its parabolic 
approximation with x = and x = — 1 are used for the uniform liquid core with 
p > Pt- Assuming that the core consists of only the npe matter without possible 
new degrees of freedom or phase transitions at high densities, the PA leads to a 
larger radius for a fixed mass compared to the full MDI EOS. Furthermore, using 
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Fig. 3. (Color online) The mass-radius relation M-R of static neutron stars from the full EOS and 
its parabolic approximation in the MDI interaction with x = and x = — 1. For the Vela pulsar, 
the const rain t of AI/I > 0.014 limits the allowed masses and radii. See text for details. Taken 
from Ref. 12ZI with small modifications. 

the full MDI EOS with x — and x = — 1 constrained by the heavy-ion reaction 
experiments, the radius of a canonical neutron star of IAMq is tightly constrained 
within 11.9 km to 13.2 km. 

5. Summery 

Using the MDI EOS of neutron-rich nuclear matter constrained by recent isospin 
diffusion data from heavy-ion reactions in the same sub-saturation density range as 
the neutron star crust, we have determined the density and pressure at the inner 
edge, that separates the liquid core from the solid crust of neutron stars, to be 0.040 
fm~ 3 < p t < 0.065 fm- 3 and 0.01 MeV/fm 3 < P t < 0.26 MeV/fm 3 , respectively. 
These constraints have allowed us to determine an improved mass-radius relation 
for neutron stars. Furthermore, we have found that the widely used parabolic ap- 
proximation to the EOS of asymmetric nuclear matter leads to significantly higher 
core-crust transition densities and pressures, especially for stiff nuclear symmetry 
energies. 
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